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Abstract 

We investigate the critical behaviour of a three-dimensional lattice yJJ 03 
model in the chiral limit. The model consists of a staggered fermion held, 
a U(l) gauge held (with coupling parameter (3) and a complex scalar held 
(with hopping parameter k). Two different methods are used: 1) fits of the 
chiral condensate and the mass of the neutral unconhned composite fermion 
to an equation of state and 2) hnite size scaling investigations of the Lee- Yang 
zeros of the partition function in the complex fermion mass plane. For strong 
gauge coupling ((3 < 1) the critical exponents for the chiral phase transition 
are determined. We find strong indications that the chiral phase transition 
is in one universality class in this (3 interval: that of the three-dimensional 
Gross-Neveu model with two fermions. Thus the continuum limit of the 
xU(p3 model dehnes here a nonperturbatively renormalizable gauge theory 
with dynamical mass generation. At weak gauge coupling and small k, we 
explore a region in which the mass in the neutral fermion channel is large 
but the chiral condensate on hnite lattices very small. If it does not vanish 
in the infinite volume limit, then a continuum limit with massive unconfined 
fermion might be possible in this region, too. 
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1 Introduction 



Strongly coupled gauge theories are interesting candidates for new mass generating 
mechanisms because they tend to break chiral symmetry dynamically. However, 
the fermions which acquire mass through this mechanism usually get confined. It 
was pointed out |l] that this can be avoided in a class of chiral symmetric strongly 
coupled gauge theories on the lattice in which the gauge charge of the fermion is 
shielded by a scalar field of the same charge. The question is, whether these models 
are nonperturbatively renormalizable at strong gauge coupling such that the lattice 
cutoff can be removed. If so, the resulting theory might be applicable in continuum 
and constitute a possible alternative to the Higgs mechanism 0. 

In this work we investigate such a lattice model in three dimensions with a vec- 
torlike U(l) gauge symmetry, which we call xU<t>3 model. It consists of a staggered 
fermion field x with a global U(l) chiral symmetry, a gauge field U G U(l) living on 
the lattice links of length a and a complex scalar field <fi with frozen length \<p\ = 1. 
It is characterized by the dimensionless gauge coupling parameter (3 (proportional 
to the inverse squared coupling constant), the hopping parameter k of the scalar 
field and the bare fermion mass am®. The unconfined fermion is the composite 
state F = <f>'x- m a phase with broken chiral symmetry, it has nonvanishing mass 
arriF in the chiral limit mo = 0. The xU<p3 model can be seen either as a general- 
ization of three-dimensional compact QED with a charged scalar field added or as 
three-dimensional U(l) Higgs model with added fermions. 

The same model has also been investigated in two and four dimensions. In 
two dimensions it seems to be in the universality class of the Gross-Neveu model 
H at least for strong gauge coupling, thus being renormalizable. Therefore the 
shielded gauge-charge mechanism of dynamical mass generation suggested in |L| 
works in two dimensions and its long range behaviour is equivalent to the four 
fermion theory. In four dimensions there is also a region in (3 (0 < (3 < 0.64) in 
which the model behaves in a very similar manner to the corresponding four-fermion 
theory, the Nambu-Jona-Lasinio model with a massive fermion whose mass scales 
at the critical point [|J. Here both models belong to the same universality class 
and have the same renormalizability properties. But for intermediate coupling 
there evidently exists a special point. It is a tricritical point at which, together 
with the composite fermion F, scaling of a particular scalar state was found. This 
composite scalar can be interpreted as a gauge ball mixing with a (j^-cf) state. Thus 
the gauge degrees of freedom play an important dynamical role and the model 
belongs to a new universality class of models with dynamical mass generation, 
whose renormalizability is of much interest j|, ||. 

In this paper we investigate the phase diagram and the critical behaviour of the 
model in three dimensions. We find that in the chiral limit mo = the xU4>3 model 
has three regions in the (3 — k plane with different properties with respect to the 
chiral symmetry. They are indicated in Fig. |l[ The region at strong gauge coupling 
(small (3) and small k is the Nambu phase where the chiral symmetry is broken 
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Figure 1: Phase diagram of the xU<p3 model for mo = 0. For (3 < 1 a clear phase 
transition between the Nambu and Higgs phases can be observed. Whether the 
region X at large (3 and small k forms a third phase or belongs to one of the other 
phases, separated from it only by a crossover, is discussed in the text. All phase 
transitions seem to be 2 nd order. 



and the neutral fermion F is massive. At large k chiral symmetry is restored and 
the fermion F is massless. This phase is labelled the Higgs phase because of its 
properties in the weak coupling limit. The third is the X region at large (3 and 
small k. It is conceivable that this region is analytically connected with either the 
Nambu or Higgs phase but it may well be a separate phase. In this region the mass 
measured in the fermion channel is large, but the chiral condensate is very small 
(within our numerical accuracy consistent with zero). 

The main result of our paper is the determination of the critical behaviour at 
strong gauge coupling. We find strong indications that the chiral phase transition 
between the Nambu and Higgs phases is in one universality class for all (3 < 1. 
It is the class of the three-dimensional Gross-Neveu model which is known to be 
(nonpertubatively) renormalizable ||. That model is the (3 = limit of the xU<t>3 
model 0. This universality means that the continuum limit of the xU<$3 model 
defines a nonperturbatively renormalisable gauge theory in which the fermion mass 
is generated dynamically by the shielded gauge-charge mechanism. However, it also 
means that in this (3 region the gauge field is auxiliary and the xU<t>3 model does 
not represent a new class of field theories. 

The chiral properties of the region X are elusive and their determination would 
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require substantial effort and resources. This is beyond the scope of the present work 
and we made only an exploratory investigation. But we point out that, provided 
the chiral symmetry is broken there, the phase transition between the region X 
and the Higgs phase gives rise to another possible construction for a continuum 
theory containing an unconfined fermion with dynamically generated mass. It could 
continue to be in the universality class of the three-dimensional Gross- Neveu model. 
But experience [|], [5| with the four-dimensional model in the vicinity of the tricritical 
point suggests that at larger j3 the gauge degrees of freedom are dynamical and a 
new universality class may be present. This interesting possibility, and the possible 
existence of a tricritical point in three dimensions, deserves further study. 

Our investigation is mainly based on two methods: firstly, via fits to an equation 
of state and, secondly, via a finite size scaling investigation of the Lee- Yang zeros 
in the complex fermion mass plane. The investigation of a phase transition via 
fits to an equation of state is quite reliable because the finite size effects we find 
close to the phase transition are usually small. Therefore we expect a simple finite 
size scaling, describe it by an empirical formula and extrapolate observables to the 
infinite volume. Then we do a simultaneous fit to the fermion mass amp and the 
chiral condensate 

As first pointed out by Lee and Yang || ||, the determination of the finite size 
scaling behaviour of the complex zeros of a partition function could be a direct 
method for the determination of the critical properties of the associated theory. 
In this paper we investigate these zeros of the canonical partition function in the 
complex bare fermion mass plane. These zeros control the fermion condensate and 
its associated susceptibilities [H], [11], physical quantities which are often measured 



directly on the lattice and used, via finite size scaling, to determine the critical 
behaviour. 

In the region X, where the chiral condensate is very small, both methods fail 
to provide reliable results. A small condensate suggests that the Lee- Yang zeros 
cannot be near to the physical region. Nevertheless, it is of interest to investigate if 
the closest zeros can be determined with sufficient accuracy to ascertain their finite 
size scaling (and hence that of the condensate). 

The paper is organized as follows. In the next section we introduce the model 
in detail, define the observables we use and briefly summarize the method of the 
Lee- Yang zeros. In section |3| we present evidence for the universality at strong gauge 
coupling. In section f| we present the results obtained at weak coupling and discuss 
their possible interpretations. In the last section our results are summarized. 

2 The model 

The xU(f>3 model is defined on a 3-dimensional cubic lattice with periodic boundary 
conditions except for antiperiodic boundary conditions for the fermion field in the 
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"time" direction. The action reads: 



with 



S x U(t> — S x + Su + S^, (1) 

1 3 

S x = g Yl X* Yl Vijai{ u x^Xx+n ~ Ul_^Xx-n) + am x x Xx , 

X A* = l 33 

= /?^(1-Re[/^), 

3 

Here Xx are the Kogut-Susskind fermion fields with t]^ x = (— l) Xl ^ Be- 
cause of doubling our model describes two four-component fermions (Nf = 2). The 
bare mass arriQ of the fermion is introduced for technical reasons. We are interested 
in the chiral limit mo = 0. The a in front of mo indicates that we have to distinguish 
between the chiral limit in the continuum (mo = 0) and the continuum limit of the 
lattice model, where am — > can also be achieved by a — > at nonzero m . 

Ux,n represents the compact link variable and U X)liV is the plaquette product of 
the link variables U Xifl . 

The hopping parameter k vanishes, if the square of the bare mass of the scalar 
field is +oo, and is infinite if the bare mass squared is — oo. The scalar field (f) 
has frozen length \<f>\ = 1. This choice is made in order to restrict the number of 
parameters of the model. Without that, symmetries and dimensionality of couplings 
would allow several other terms in the action. 

We stress that the charges of the fundamental fields exclude a direct Yukawa 
coupling between the fundamental fields. 

The model has some interesting limiting cases. For k = the scalar field decou- 
ples and the model is equivalent to three-dimensional compact QED with fermions. 



It is known |L2| [13|, |14| [T5[ that pure compact QED has no phase transition and, as 
(3 — > oo, it is confining via a linear potential with an exponentially decreasing string 
tension. There is an indication that, with fermions, chiral symmetry is broken at 



large coupling, but at weak coupling results are inconclusive pi. It has been sug- 



gested that, in noncompact QED with fermions, the phase diagram is dependent 
on the number of flavors and that, at weak coupling, chiral symmetry is broken 



only for a small number (less than about 3.5) of fermions [|17], 18|, [L9|, ^0|. (A recent 
description of the status of three-dimensional QED can be found in ^2|.) It is 
quite probable that, at weak coupling, both the compact and non-compact formu- 
lations have quite similar properties. If so, then these (uncertain) features suggest 
that the chiral symmetry is broken in the k = limit of the phase X and thus 
presumably in the whole phase X. 
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In the weak gauge coupling limit, (3 = oo, the fermions are free with mass awio, 
and reduces to the XY 3 model. It has a phase transition at k ~ 0.27. 

At amo = oo the model reduces to the three-dimensional compact U(l) Higgs 



model. For its recent investigation with numerous references see ||23 . 

For (3 = the gauge and scalar fields can be integrated out exactly |7j and one 
ends up with a lattice version of the three-dimensional four-fermion model 

S 4f = ~Yl ^2i G XxXxX x+ ^Xx +l , - 7>V^(XxXx+n - Xx+^Xx)] + XxXx, (2) 

X ^t=l X 

the parameters G and r being related to k [7|, |T[| . 

We refer to this model as the Gross-Neveu model. Some caution is in place, 
however. There is some uncertainty whether the four-fermion model @) is a lattice 
version of the Gross-Neveu model or of the Thirring model. The four-fermion 
action @ was used in four dimensions for the study of the Nambu-Jona-Lasinio 
model e.g. in [^4, ^5], which would correspond to the Gross-Neveu model in three 



dimensions. Recently in [^, ^7] the four-fermion action (0) in three dimensions 



is interpreted as the Thirring model and similar interpretation is implied by K.- 
I. Kondo p8f . For our number of fermions, Nf = 2, the distinction might be 
unimportant and both models might actually coincide^. The Gross-Neveu model 
has a chiral phase transition and is nonperturbatively renormalisable (see || 29 



and references therein). The properties of the Nf = 2 Thirring model appear to 
be similar |2^, [26], ^7|. For our purposes the important property of the three- 
dimensional four-fermion model obtained in the (3 = limit of the xU<p3 model is 
its renormalisability, which presumably holds for both interpretations. 



2.1 Observables 

Because we are interested in the chiral properties of the model we concentrate on 
the chiral condensate and the fermion mass. 
The chiral condensate is defined by 

(xx) = <TrM" 1 > (3) 

where M is the fermion matrix. The trace is measured with a gaussian estimator. 

The physical fermion of the xU<p3 model is the gauge invariant composite fermion 
F = <$X- We measure its mass artiF in momentum space with the usual procedure, 
as described (for the model in 4 dimensions) in 0. We checked that the results are in 
good agreement with the fits done in configuration space. In the three-dimensional 
model we find the fits to G 1 to be the most stable, so we use them for the data 
shown in this paper. 

1 J.J. thanks M. Gockeler, S.J. Hands, and K.-I. Kondo for discussions on these questions. Some 
of them are exposed in pTj 
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Both observables need to be extrapolated to infinite volume. This procedure is 
described in section |3~TT| . 

We remark that the required numerical effort for the study of the xU4>3 model 
was very high. We needed significantly more matrix inversions than for the four- 
dimensional case || at the same volume and amo. Their number also depends 
significantly on (3: The simulation at (3 — required about 0(1000) conjugate- 
gradient steps, about 2-5 times more than at larger /3-values. Surprisingly, the 
number of required steps scattered in a very broad interval. The maximal step- 
number was at least a factor of 2-3 above the average. This might be connected 
with the observation, that the chiral condensate has a very asymmetric distribution. 

2.2 Equation of State 

A standard way to analyse the critical exponents of a chiral phase transition is via 
the use of an equation of state (EOS). Normally data close to the phase transition 
can be well described by such an ansatz. In our model for fixed f3 this equation 
reads 

am = (xx) 5 F ((k - k c ){xx)' 1,Px ) , F(x) = Rx + S . (4) 

Here (xx) is the infinite volume value of the chiral condensate for given amo, K an d 
j3. f3 x and 8 are the exponents defined in analogy to a magnetic transition. The 
index x is added to distinguish the exponent and the coupling. The scaling function 
F is used in its linear approximation and R and S are free constants. We apply 
this equation in the region for which k ~ k c and where we might expect the scaling 
deviations to be small. 

It is also useful to assume the corresponding scaling equation for the fermion 
mass: 

am = (am F ) 1/9 G ((« - K c )(am F y l,v ) , G(x) = Ax + B . (5) 

The exponent v is the correlation length critical exponent in the chiral plane (am = 
0). v is an analogous exponent obtained if one approaches the critical point from 
outside the chiral plane. The two exponents have to be distinguished. At fixed f3 in 
the chiral plane (amo — 0) the fermion mass scales with u: amp oc (k — n c Y\am =o, 
whereas for all other straight paths into the critical point (for example k = k c ) it 
scales with am as: amp oc am u \ K - Kc ocarn - This is indicated in figure H in the plane 
(3 = const. 

Figure ^| also illustrates that for k < k c the chiral condensate changes sign and 
makes a jump if one crosses the line amo = 0. This means that it is a line of first 
order phase transitions. For k > k c the line becomes a line of second order phase 
transitions on which the fermion mass gets critical. In between there is a critical 
point (k = k c ). 

If hyperscaling holds, only two of the four exponents defined by the equations 
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Figure 2: First order chiral phase transition line (bold line) and the critical point in 
a plane (5 = const. The fermion mass amp scales with exponent v in the direction 
tangential to the transition line and with exponent v in any other direction. 



of state are independent. The corresponding scaling relations are: 

1 / 1\ 

5 = and B v — v d , 

dv-\ x \ v ) 

where d — 3 is the space-time dimension. 



(6) 



2.3 Lee- Yang zeros 

The canonical partition function, after integration over the Grassmann variables 
and using the irrelevance of overall multiplicative factors, can be defined as: 



Z(/3, k, amo) 



J dUd(f> det M[am , U]e 
J dU dtp det M[am ,C/]e' 



U4> 



(7) 



Here Su<j> = Su + S<p, M is the usual fermionic matrix for Kogut-Susskind fermions 
and om is some (arbitrary) "updating" fermion mass at which the ensemble of 
gauge fields is generated. 

Since the mass dependence of M is purely diagonal, the partition function can 
be written as the average over the ensemble of the characteristic polynomials of M, 
i. e.: 



Z(/3, K,am ) = 



^=oCn[U[(3,K}}(am ) 2 n 
det M[am , U[/3, k}} 



n=0 



amo 



(8) 
(9) 
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The coefficients C n of the characteristic polynomial are obtained from the eigen- 
values of M[0, U] which are imaginary and appear in complex conjugate pairs. In 
the simulations described below they were obtained using the Lanczos algorithm. 

The Lee- Yang zeros are the zeros of this polynomial representation of the parti- 
tion function. The zeros were found by using a standard root finding algorithm on 
the equivalent sets of polynomials generated as in Eq. (8): 

v 

Y^^aml-amiY (10) 

n=0 

for a set of am, in the region where we expect the lowest zeros to occur. This 
allowed us to avoid the problems associated with rounding errors in the root-finder. 
We required that a zero be found consistently for the subset of the arhi closest to 
it. These zeros in the bare mass we label as yi in the following. 

The errors in the Lee- Yang zeros are estimated by a Jacknife method. The 
coefficents for each lattice size were averaged to produce 6 subsets of averaged 
coefficients, each taking into account 5/6 of the measurements. These 6 different 
sets of coefficients give 6 different results for the Lee- Yang zeros from which the 
variance was calculated. 

The critical properties of the system are determined by the zeros lying closest 
to the real axis. The zero with the smallest imaginary part we label y\. It is also 
called edge singularity. With increasing finite volume it converges to the critical 
point. For a continuous phase transition the position of the zeros closest to the real 
axis in the complex plane is ruled by the scaling law 

Vi {p, K, L) - y R (fi, k, oo) = AiL~ 1/s , (11) 

where the A^s are complex numbers. The exponent s = s(/3, k) describes the finite 
size scaling of the correlation length. For our model Ur((3, k, oo) = and we ignore 
it in the following. 

It immediately follows that the real and the imaginary parts of the zeros should 
scale independently with the same exponent. In particular, for the zero y\ closest 
to the chiral phase transition (at am Q = 0) 

Imy 1 (f3,K,L) = A I L- 1 / s , (12) 

with a similar scaling behaviour for Rej/i(/3, k, L) via Ar. In practice the real part 
of the zero is much smaller than its imaginary part or is identically zero. So eq. fll2]) 
usually provides a more reliable measure of the exponent than the scaling of the 
real part. 

Although the above scaling law was originally established for the case of a con- 
tinuous phase transition, it can also be extended to that of a first order phase tran- 
sition. Since there is no divergent correlation length, the exponent is determined 
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Figure 3: Schematic plot of the exponent s in infinite volume (full lines and dot) 
and the effective s in finite volume (dashed line). 
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Figure 4: Expected finite size scaling of the zero y\ with lattice size. 



only by the actual dimension of the system. In this case, for a three-dimensional 
model we expect s = |. 

At the critical point (k = k c ) we expect s to be equal to u, because the fermion 
correlation length should be the relevant one. In the symmetric phase (k > k c ) 
we expect scaling with s — 1, because amp oc arriQ. This behaviour is indicated in 
Fig. |H by the full lines and the dot. 

In practice it is important to understand the scaling deviations on a finite lattice. 
The expected behaviour is shown schematically in Fig. [|. Far away from the critical 
point we expect linear scaling in the log- log plot with s = 1/3 in the broken phase, 
and s — 1 in the symmetric phase. At the critical point we expect linear scaling 
and the exponent should be s = v. These expectations are indicated by the full 
lines. Close to the phase transition, we expect a crossover. For small lattice sizes 
the exponent should be close to v and then change to 1/3 and 1, respectively, if the 
lattice size is increased and the true scaling shows up. This is indicated in Fig. |] by 
the dashed lines. For a set of lattice sizes this defines an effective s which smoothly 
goes through v at the critical point. Such an effective s is represented in Fig. [I] by 
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a dashed line. 

Therefore, in order that the critical exponent can be determined, we must either 
know the position of the critical point accurately or have many simulations on large 
lattices so that the scaling deviations can be measured accurately. In practice the 
limited knowledge of the position of the critical point leads to the largest uncertainty 
in the determination of v by this method. 



3 Universality at Strong Coupling 

At strong coupling the chiral phase transition can be seen clearly and we investigate 
the scaling behaviour and the universality along this line. 

We determined the Lee- Yang zeros, the chiral condensate and the fermion mass 
for various values of k at /3 = 0.00 and 0.80. In this section we want to investigate 
how the transition changes as (3 increases from zero. Therefore we have investigated 
the scaling of the data at (3 = 0, i.e. the four-fermion model, as a reference and 
compare it with the scaling found at (3 = 0.80. 



3.1 Equation of State 

Here we determine the critical exponents of the chiral phase transition by using the 
EOS for arriF and (xx}. Although we did simulations on lattices up to 24 3 , our 
conclusions still depend to some extent on our choice of ansatz for the extrapolation 
of arriF and (xx) t° infinite volume. Fig. |5| shows, as an example, our data for amp 
at (3 = 0.80 and am = 0.01 plotted against 1/L 2 . 
For the extrapolation we tried three approaches: 

ampiL) = amp( oo) + A— , (13) 
amp(L) = amp(oo) + A— , (14) 
am F (L) = am F (oo) + A— exp(— amp(oo) L) . (15) 

Each has two free parameters: am^ (oo) and A. To judge the quality of the fits 
we first compared the x 2 P er degree of freedom using our data on 16 3 , 20 3 and 24 3 
lattices. This was done at the values of (3 and k at which we have good statistics. 
Our results are shown in table [I]. It turned out that the results for these lattice 
sizes are not conclusive as to which extrapolation formula should be used, because, 
for each ansatz, all x 2 P er degree of freedom are usually below 1. 

However, the fit with eq. (|13|) is significantly prefered if compared with the 12 3 
data. We therefore adopted this fit for our extrapolations, but data from the 12 3 
lattice was not included. Such an extrapolation is indicated in Fig. || by the dotted 
lines. 
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Figure 5: Data for amp at (3 = 0.80 and arriQ = 0.01 plotted against 1/L 2 . The 
dotted lines are a fit with eq. (pT3|) to the data with L > 16. 



/? k amo 


Fit 1 

amp x 2 


Fit 2 

arriF \ 2 


Fit 3 

arriF \ 2 


0.00 0.95 0.01 
0.00 1.00 0.01 
0.00 1.05 0.01 
0.00 0.95 0.02 


0.269(3) 0.30 
0.194(4) 0.22 
0.132(5) 0.23 
0.337(3) 0.41 


0.256(7) 0.53 
0.180(9) 0.45 
0.111(10) 0.04 
0.327(6) 0.22 


0.278(2) 0.01 
0.202(2) 0.01 
0.142(3) 0.62 
0.343(1) 2.07 


0.80 0.40 0.01 
0.80 0.42 0.01 
0.80 0.43 0.01 
0.80 0.45 0.01 


0.441(8) 0.12 
0.276(8) 0.01 
0.215(5) 0.07 
0.122(2) 0.01 


0.437(2) 0.15 
0.247(16) 0.08 
0.191(9) 0.01 
0.087(1) 0.14 


0.445(3) 0.03 
0.295(4) 0.31 
0.230(2) 1.33 
0.137(3) 0.36 



Table 1: Results of the fits to the finite size behaviour at different couplings and 
masses on 16 3 , 20 3 , 24 3 lattices. The extrapolated infinite volume mass arriF = 
amp (oo) and the y 2 per degree of freedom for the three fits are given: with eq. fliTf ) 
(Fit 1), eq. ([H (Fit 2) and eq. (Fit 3). 
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a 




V 


V 


A 


B x 2 


amp: 


0.00 
0.80 


0.987(33) 
0.425(5) 


0.91(22) 
0.78(14) 


0.43(8) 
0.40(5) 


1.1(3) 
3.6(6) 


0.38(8) 0.72 
0.33(5) 0.82 







K c 


Px 


5 


R 


s X 2 


(xx)- 


0.00 
0.80 


0.983(12) 
0.429(7) 


0.56(5) 
0.56(10) 


3.1(3) 
3.0(5) 


1.3(2) 
4.4(12) 


1.1(3) 0.84 
2.4(13) 0.58 



Table 2: Results of fits at (3 = 0.00 and 0.80 using the equations of state. The 
upper table shows the results of the fit of am F based on eq. (|5|), the lower table 
those of (xx) based on eq. (f|). 



For the chiral condensate the finite size effects are in general smaller and with 
opposite sign. Again, consideration of the 12 3 lattices favoured a fit ansatz analo- 
gous to eq. (|T3|). 

We describe in detail the analysis in which we used the ansatz of eq. fll3|) to 
extrapolate all our data for arriF and (xx)i obtained on 16 3 and larger lattices, to in- 
finite volume. The error was calculated with the MINOS routine from the MINUIT 
library. All results presented in the following change somewhat quantitatively, but 
not qualitatively, if a different extrapolation formula is used. 

The data at different k and amo, extrapolated to the infinite volume, were 
analyzed by means of the EOS. We included only the data at am = 0.01 and 
arriQ = 0.02. The chosen k range was 0.80 . . . 1.05 for (3 = 0.00 and 0.38 . . . 0.47 for 
(3 = 0.80. 

As a first step we analyzed the data for amF and (xx) independently and fitted 
to their corresponding EOS ([5]) and (Q). The results are given in table 0. As can 
be seen, for both /3's the critical k values Kjq £1X6 identical within the error bars. 

As a next step we performed a simultaneous fit with one common k c for amF 
and (xx) (table ^j. A very good fit to all the data was obtained. 

Then we checked the scaling relations (0). Calculating (3 and 8 with v and v 
gives (3 X = 0.54(20) and 5 = 3.8(15) for (3 = 0.00 and (3 X = 0.39(25) and 5 = 5(3) 
for j3 = 0.80. The agreement with the fit is quite good. Note that in (^), dv = 3u 
is close to 1 and hence the statistical errors are increased. 

We also tried a third fit in which we assumed the validity of the scaling relations 
(0). The result is shown in Figs. |6] and |7| and summarized in table (|. As one can see, 
the quality of the fit is still good and x 2 are reasonable. The figures also show the 
prediction of our fit for the fermion mass and chiral condensate at amo = 0.04 and 
0.06. Only small deviations are visible. We therefore conclude that @ is consistent 
with our data. 

The values of the exponents z/, z>, and (3 X in table |2] and table § agree with those 
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Figure 6: (a) Fermion mass and (b) chiral condensate for f3 = 0.00. The data 
are our extrapolation into the infinite volume. The fit assumes the validity of the 
scaling relations and is described in the text. The parameters are given in table f|. 
The dashed line shows the extrapolation into the chiral limit. 
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p 




V 


V 


Px 


5 


x 2 


0.00 
0.80 


0.983(12) 
0.425(4) 


0.88(8) 
0.78(11) 


0.42(3) 
0.40(4) 


0.56(5) 
0.47(5) 


3.1(3) 
3.4(3) 


0.71 
0.71 



Table 3: Results of fits of amp and (xx) at /5 = 0.00 and 0.80, using both equations 
of state with a common k c . 



P 


k c v v x 2 


Px S 


0.00 
0.80 


0.981(6) 0.79(2) 0.437(5) 2.2 
0.425(2) 0.75(2) 0.431(6) 2.3 


0.56(4) 3.2(2) 
0.51(4) 3.4(2) 



Table 4: Results of our fits using the equations of state at (3 = 0.00 and 0.80 with 
one k c and the scaling relations @ at (3 = 0.00 and 0.80. 



in table |]. Thus all three fitting procedures gave consistent results at each (3. 

Furthermore, the exponents obtained at (3 = 0.00 and f3 = 0.80 agree within 
errors. This is a strong signal that the chiral phase transition is in one universality 
class at these /3's and probably also for those in between. The difference in the 
exponents of the last fit, which is somewhat larger than the pure statistical errors 
may for example be the result of small scaling deviations. 

To estimate the uncertainty due to the choice of the extrapolation formula, we 
repeated the above procedures using the extrapolation fll!]) . The results are given 
in table || The x 2 's are even smaller and the exponents differ by a little more than 
one standard deviation. Although the agreement for the two /3's is less good, it is 
still compatible with universality if one takes into account that the error bars only 
reflect the statistical errors and not the uncertainty due to scaling deviations. 



p 


K c V V X 2 


Px 6 


0.00 
0.80 


0.968(9) 0.76(3) 0.426(8) 0.89 
0.419(3) 0.66(4) 0.409(10) 0.64 


0.50(6) 3.6(3) 
0.37(6) 4.4(6) 



Table 5: Results of fits, when we use eq. ( |14"D for the extrapolation to infinite volume. 
As in table |], the equations of state with common n c and the scaling relations (|(J) 
are used. 
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Figure 8: Imaginary part of zero y\ as function of the lattice size for different k at 
(3 = 0.00. The different straight lines should help to investigate the linearity and 
claim only for k = 0.975 ~ k c (full line) to describe the data well. 



3.2 The Finite Size Scaling of the Lee- Yang Zeros 

The Lee- Yang zeros were found to be purely imaginary at all values of k in the 
strong coupling region. For small k they are equally spaced, consistent with a 
strong first order transition in the condensate as the bare mass amo goes through 
zero. 

Figs. |8| and |9| show the finite size scaling behaviour of the edge singularity at 
various k for (3 = 0.00 and 0.80, respectively. Our data confirm the expectations 
presented in section O] and Fig. 01. Close to the critical point we see the expected 



crossover: for small lattices the exponent is close to v and shifts for increasing 
lattice size to the exponent 1/3 or 1. 

At small k the exponent s is consistent with a first order phase transition, 
s > |. At large k is s < 1. Close to the critical point, determined in the previous 
section and given in table [|, the data scale linearly on the log-log plot allowing 
determination of v. 

At the k points closest to k c we expect s ~ v. At (3 = 0.00 we find at re = 
0.975 ~ k c = 0.981(6) the exponent s = 0.440(4) in excellent agreement with 
v = 0.437(5), as determined in the previous section. For (3 = 0.80 we did the 
simulations at k = 0.43 slightly larger than the k c = 0.425(2) obtained from the 
EOS. Not unexpectedly we found s = 0.447(5), slightly larger than v = 0.431(6) 
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,3 = 0.80 

fit: A*L~ 1/S : 

*:/c = 0.70 
s= 0.865(17) 

t:/c = 0.60 
s= 0.853(33) 

Q:»c = 0.50 
s= 0.593(11) 

■»:k = 0.47 
s= 0.5396(80) 

□ :k = 0.45 
s= 0.5007(58) 

&:k = 0A3 
3= 0.4468(43) 

*:(C = 0.40 
s= 0.374(21) 

0:^ = 0.30 
s= 0.3417(29) 

Figure 9: Imaginary part of zero y\ as function of the lattice size for different k at 
= 0.80. 



from the EOS. This shows the great importance of the knowledge of the critical point 
for a precision measurement of v. Within these uncertainties the agreement is very 
good and again confirms the independence of v from f3 and hence the universality. 

If one analyses these plots without the knowledge of the critical point determined 
with the EOS, the critical point can also be determined by looking for linearity of 
ln[Im(yi)] as a function of InL. For this purpose we show in Figs. |10| and [ll| the 
quantity ln[Jm(?/i)] + l/i>lnL. The addition of the second term makes the plots 
approximately horizontal and so allows us to enhance the vertical scale making the 
error bars clearer. In the figures we have used v = 0.440 and 0.501, respectively. 
The dashed lines are a linear extrapolation of the data points at the two lowest 
L-values. They provide a guide as to the linearity of the data. 

Figs. [XOj and [U] suggest a larger k c and s(k c ) than those obtained from the EOS 
analysis. For example, at (3 = 0.80, Fig. |TT] would suggest k = 0.45 as the point 
closest to k c , with v ~ s(0.45) = 0.501(6). However, the difference between the two 
methods of analysis is about 10% which is the same size as the statistical error. 

All in all, this demonstrates the reliability of the methods we have used. Both 
(very different) methods agree rather well and their combination is very useful. 
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3.3 Universality at strong coupling 



Our data are a good indication that the chiral phase transition of the xU(f>3 model is 
in the same universality class at (3 = 0.00 and (3 = 0.80. Assuming this universality 
we combine the results for exponents at both (3 values and determine the exponents 
of this chiral phase transition to be v = 0.75(10) and v = 0.43(2). The errors take 
into account the uncertainties discussed above. These values of v and v correspond 
to (3 X = 0.51(11) and 5 = 3.45(71). We note that the position of the critical point 
at (3 = 0, as well as the results for S and j3 are compatible^ with those obtained 
for the N f = 2 case in J27| (fi x = 0.57(2) and 5 = 2.75(9)). In that work the same 
action fl2|) has been simulated, though in a somewhat different representation by 
means of auxiliary fields than the (3 = limit of the xUcj)^ model. 

It is very likely that the chiral phase transition is in the same universality class 
for (3 between and the onset of the X region around /3 ~ 1. The universality might 
be expected at small (3 because of the convergence of the strong coupling expansion. 
But our data are (to our knowledge) the first indication that this is true for a large 
f3 interval. 

This result indicates that the xU<p3 model is renormalizable in this region of 
f3. So it is a nontrivial example in three dimensions for the shielded gauge-charge 
mechanism of fermion mass generation proposed in JT|. 

We note that, in the scaling investigation, the chiral condensate, which is a pure 
fermionic operator, and the mass of the fermion F, which is a combination of the 
fermion and the scalar field, have been used. Both seem to scale in a way which 
can be well described by the usual scaling relations. 

The universality on the other hand also means that, with respect to the three- 
dimensional Gross-Neveu model, nothing substantially new happens at small and 
intermediate (3 and no new physics arises on scales much below the cutoff. The 
scalar field shields the fermion x giving rise to the fermion F equivalent to the 
fermion of the four-fermion theory. We find no indication that composite states 
consisting only of fundamental scalars or gauge fields, which would not fit into the 
Gross-Neveu model, scale at the chiral phase transition. 

The bosonic fields appear to be auxiliary at strong gauge coupling, as they are 
in a rigorous sense at (3 — 0. As indicated by the results in four dimensions 
H, §], this may change as the gauge coupling gets weaker. We therefore performed 
some studies at larger values of (3. The results are described in the next section. 

4 Explorative study of the weak coupling region 



2 We thank S.J. Hands for pointing out to us this compatibility. 
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Figure 12: (a) fermion mass and (b) chiral condensate for different small arriQ as 
function of (3 at k = 0.25 on the 8 2 16 lattice. 



4.1 Condensate and fermion mass 

We should point out again that the compact three-dimensional QED with fermions 
is not fully understood at large (3. It is not clear, at large (3, if there is chiral 
symmetry breaking and confinement via a linear potential. This uncertainty extends 
also to our model with k small. A clarification of these difficult questions would 
require a substantial effort far beyond the scope of this paper. Thus our aim is to 
perform an explorative study only and to get some insight into this as yet unexplored 
region. Also we want to see how far the methods applied successfully at strong 
coupling can be of use also at weaker coupling. The physical interpretation of our 
results will leave room for several scenarios. 

Fig. [12] shows the fermion mass and condensate at k — 0.25 as a function of j3, 
at three values of the bare fermion mass. The neutral fermion mass decreases for 
increaing f3 but then stabilizes with amp > 1. So it is clearly nonzero at all f3 and 
again only weakly dependent on the bare fermion mass. The condensate is large at 
small (3 (the Nambu phase) but rapidly decreases with (3 and becomes very small 
(zero?) in the chiral limit for (3 > f3x — 1-3. Thus here a new, weakly coupled 
region is encountered. 

In order to see how this region is related to the Higgs phase at large k, the 
fermion massQ and condensate are shown in Fig. 13 at (3 = 2.0 and three values of 



3 We shold remark that for large (3 the agreement of the different fits for the fermion mass is 
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Figure 13: (a) fermion mass and (b) chiral condensate for different small amo as 
function of k at (3 = 2.00 on the 8 2 16 lattice. 



the bare fermion mass. For nonvanishing bare mass, where the simulations have 
been performed, the mass of the neutral fermion is large for k < Kx — 0.27 whereas 
it is small for larger k. It is only very weakly dependent on the bare fermion mass 
and therefore we expect this behaviour to persist in the chiral limit. For k > Kx its 
small nonzero value probably vanishes in the infinite lattice size limit. 

The condensate, as expected, does depend strongly on the bare mass amo but 
does show a crossover behaviour at n = k x with {xx) ( am o) k<k x < (xx)( am o)«>Kx- 
However, at large k we believe that we are in the Higgs phase where the condensate 
is zero in the chiral limit. It is therefore conceivable that, in this limit, it is zero at 
f3 = 2.0 for all k. It is very surprising, however, that, at fixed bare fermion mass 
and lattice size, the condensate tends to slightly increase with increasing k, quite 
in contrast from its behaviour in the strong coupling region. Of course, this can 
change in the infinite volume and chiral limit. 

At the coupling k = k^, where the fermion mass shows a possible crossover 
behaviour, there is also a peak in the susceptibility of the link energy which increases 
with increasing lattice volume. However, a careful finite size scaling analysis would 
be needed to determine if it indicates a phase transition or a crossover. 



Fig. [14] a confirms the weak dependence of the fermion mass amp on the bare 



not as good as that at small /3. The qualitative behaviour is not influenced by this. This same 
feature was also observed in the four dimensional model for large f3 but is not understood up to 
now. 
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Figure 14: (a) Fermion mass and (b) chiral condensate for different k as function 
of am Q at fi = 2.00 on the 8 2 16 lattice. 



mass. At k — 0.22 below Kx, the fermion mass is large and stays clearly nonzero 
at am . Above Kx, at k — 0.34, the fermion mass is too small for the lattice size 
used and might vanish in the infinite volume limit. 

Fig. |T4| b shows that, at fi — 2.0 and k just below or above Kx, the condensate 
extrapolates linearly in omo to a very small value or zero. As we shall see below, 
this is due to the Lee- Yang edge singularity in this region being relatively distant 
from the real axis. 

A naive extrapolation to the chiral limit would thus classify the region at small 
k (k < Kx) and large (3 {(3 > fix) as a phase with zero chiral condensate and 
nonvanishing fermion mass. But the condensate could also remain very small but 
nonvanishing. Because of this uncertainty we label this region X. Its boundaries f3x 
and Kx may slightly depend on k and (3, respectively. 

It would be surprising if in the region X the fermion mass was different from 
zero with unbroken chiral symmetry. There are essentially two scenarios avoiding 
such a paradox. 

1) Chiral symmetry breaking persists at small k also for (3 > fix- X is light, 
because the chiral condensate is very small, though nonzero. F is a bound state of 
and x- The binding might be quite loose, presumably by a weak linear confining 
potential, which one expects in pure U(l) in 3d [14], [T5| . F is heavy essentially 
because <fi is heavy. The transition at (3 — fix is probably a cross-over, but a genuine 
phase transition is not excluded. 
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In this scenario the region X must be separated by a chiral phase transition at 
Kx(/3) from the Higgs phase. As our data around Kx do not indicate any metastabil- 
ity, it would be a higher order transition and a continuum limit should be possible. 
Thus an interesting continuum limit with massive unconfined fermion might exist. 

2) Chiral symmetry is restored at (3 = fix and the chiral condensate thus vanishes 
identically and x is niassless in X. The F channel gets contribution from the two- 
particle state and x- This contribution appears as a massive state because 4> is 
heavy. This state presumably cannot be a bound state in the chiral limit because 



of the old argument of Banks and Casher [|30| : fermion on a closed orbit must be 
able to flip its helicity, i.e. existence of the bound state implies chiral symmetry 
breaking. We cannot distinguish between a bound state F and a two-particle state 
<p + x looking at the F channel only (as we did). X could be connected to the Higgs 
phase, where we expect the same spectrum. 

Which of these scenarios is true might be investigated in the limit case k = 0, 
i.e. in the three-dimensional compact QED. The results in the noncompact case 



U\, [L9|, g(J m ight be applicable at weak coupling also to the compact one. As 
the number of fermions in our case is below the critical number ~ 3.5 of fermions 
in the noncompact model, the more interesting scenario 1) seems to be preferred. 

4.2 The Lee- Yang zeros at Weak Coupling 

In an attempt to clarify the situation at weak coupling we have also investigated 
the Lee- Yang zeros in the region X. The edge singularity^ y\ has a nonzero real 



part in this region. The finite size scaling of the lowest zeros is shown in Fig. 13 
for (3 = 2.00, k = 0.15, a point in the middle of the region X. The real part of the 
low lying zeros is clearly nonvanishing. Then the first two zeros have within the 
numerical precision identical imaginary part but their real parts differ by a factor 
of about 3.5. 

The first two zeros have imaginary parts so close as to be indistinguishable 
within statistical error. We have assumed continuity in the behaviour of their real 



parts as a function of lattice size when plotting Fig. |T5 . 

These imaginary parts scale linearly in the log-log plot with an exponent s = 
0.73(4) where the error is given by the difference between the two zeros. Their real 
part has somewhat larger errors but scales within the numerical precision with the 
same exponent. This pattern for the edge singularity is found throughout the region 
X. Note that this behaviour is only observed in the edge singularity. 

Fig. [TBI shows the behaviour of the imaginary part of the edge as a function 



of lattice size at k = for various (3. There is a crossover between (3 = 1.25 and 
(3 = 1.88, i.e. from the Nambu phase, where the imaginary part of the edge zero is 
small and the transition first order, to a region where the imaginary part is large 



4 The zeros must appear in conjugate pairs. We define the edge singularity in this region to be 
the zero with smallest positive imaginary part and count each pair only once. 
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Figure 15: (a) Real and (b) imaginary part of the first two Lee Yang zeros (sorted 
by their positive imaginary part) for (3 = 2.00, k — 0.15 as function of lattice size. 
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(with nonzero real part) again consistent with a vanishing condensate. No scaling 
deviations can be observed for f3 > 1.88. In the region X the critical exponent has 
a very weak dependence on (5 and increases only very slowly on further increase of 

K. 

If the scaling in the region X is different from that in the other regions, the most 
naive expectation would be that the exponent s is universal. This is compatible 
with our data at (5 ~ 2 but not at (5 — 5. Further simulations on larger lattice are 
necessary to confirm this difference. 

It is to our knowledge the first model in which scaling of the real part of the 
edge singularity to zero has been observed. We do not understand the implications 
of this behaviour. It may well be a key point in understanding the critical nature 
of region X. 

Summarizing, the region X can be distinguished from the Nambu and Higgs 
phase by the edge singularity having a real part (on a finite lattice) and a scaling 
which cannot be described by either an exponent s = 1/3 or s = 1. If the exponent s 
is different from those in Nambu and Higgs phases then the region X is presumably 
a new phase. We have been unable to determine if the chiral symmetry is broken 
or not in this region. 

5 Conclusions 

We have presented an extensive analysis of the phase structure of the three-dimensional 
fermion-gauge- scalar model. The analysis has been made possible by the applica- 
tion of two different methods: l)fits to an equation of state of the chiral condensate 
and the mass of the physical neutral fermion and 2)finite size scaling investigations 
of the Lee- Yang zeros of the partition function in the complex fermion mass plane. 

Our investigations showed that there are three regions in the f3-K plane with 
possibly different critical behaviour in the chiral limit: 

a) The region at small k and strong gauge coupling, where the chiral symmetry 
is broken and the neutral physical fermion is massive, called the Nambu phase. 

b) The region at large k, where the chiral symmetry is restored and the physical 
fermion is massless, called the Higgs phase. 

c) A third region at weak coupling and small k, where the chiral condensate is 
zero within our numerical accuracy but the neutral fermion mass is large, called the 
X region. This region can analytically be connected with either the Nambu or Higgs 
phase but it may well be a separate phase. If chiral symmetry is not broken in this 
region, then the mass observed in the fermion channel is presumably the energy 
of a two-particle state. Otherwise this region might be an interesting example of 
dynamical mass generation of unconfined fermions. If the continuum limit is taken 
at the Higgs phase transition, the gauge fields should play an important dynamical 
role and the model would not fall into the universality class of the three-dimensional 
Gross-Neveu model. A further investigation of this possibility is highly desirable. 
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At strong gauge coupling, the chiral phase transition can be clearly localized, and 
there are strong indications that it is in one universality class for all (3 < 1: that of 
the three-dimensional Gross-Neveu model, which is known to be non-pertubatively 
renormalizable. This demonstrates that the three-dimensional lattice xU4>s model 
is a nonperturbatively renormalizable quantum field theory and the shielded gauge- 
charge mechanism of fermion mass generation |I[] works in three dimensions. 
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